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Spin echo techniques are essential for achieving long coherence times in solid state quantum memories for 
light because of inhomogeneous broadening of the spin transitions. It has been suggested that unrealistic levels 
of precision for the radio frequency control pulses would be necessary for successful decoherence control at 
the quantum level. Here we study the effects of pulse imperfections in detail, using both a semi-classical and a 
fully quantum-mechanical approach. Our results show that high efficiencies and low noise-to-signal ratios can 
be achieved for the quantum memories in the single-photon regime for realistic levels of control pulse precision. 
We also analyze errors due to imperfect initial state preparation (optical pumping), showing that they are likely 
to be more important than control pulse errors in many practical circumstances. These results are crucial for 
future developments of solid state quantum memories. 

PACS numbers: 



I. INTRODUCTION 

Quantum memories for light |Q1 0] are key elements of 
quantum repeaters JH0], which are necessary to distribute en- 
tanglement over long distances for future quantum networks 
(3ft . Quantum memories based on atomic ensembles fl-SEI- 
Il4f1 are particularly attractive in practice because the light- 
matter coupling is enhanced by the large number of atoms and 
by collective interference effects. In the retrieval process col- 
lective interference can strongly enhance the re-emission of 
the stored light in a well-defined direction, compared to the 
non-directional background emission. This makes it possi- 
ble to achieve high retrieval efficiencies 

and small 

noise-to-signal ratios. 

For long-distance applications such as quantum repeaters 
it is essential for the memories to allow long storage times. 
This can be achieved by usin g lo w-l y ing atomic states (spin 
states) for storage Bll43l la-lia, Il3l Il4ll . However, spin states 
are typically affected by inhomogeneous broadening, i.e. dif- 
ferent atoms in the ensemble have slightly different energies. 
For atomic gases this can be due to residual external magnetic 
fields or intensity-dependent light shifts (for optical dipole 
traps). In atomic gases it is possible to work with field- 
insensitive clock transitions ifTa, [l9ll to suppress inhomoge- 
neous broadening due to magnetic fields. 

Solid-state atomic ensembles, such as rare-earth ion doped 
crystals, are attractive because there are no unwanted effects 
due to atomic motion and because solid-state systems promise 
enhanced scalability. However, they also have inhomoge- 
neous broadening of the spin transitions. For example, in 
rare-earth doped crystals the rare-earth ions themselves pro- 
duce a spatially varying potential due to spin-spin interactions 



Inhomogeneous broadening is important because in the ab- 
sence of control techniques it limits the coherence time of 
collective memory excitations to the inverse of the inhomo- 
geneous linewidth, which is typically in the tens of microsec- 



onds range and is much shorter than the desired storage times. 
This effect can be compensated using spin-echo techniques, 
such as the application of a single or a pair of n pulses. The co- 
herence time can be further extended even beyond the single- 
atom T2 time by applying chains of n pulses (bang-bang con- 
trol) HQ. 

In practice the control pulses are never perfect. In a recent 
experiment i23ll the most important imperfection was shown 
to be an inhomogeneity in the rf intensity across the sample, 
leading to a variation of about 1 % in the total pulse area seen 
by individual atoms. 

In Ref. i24ll the authors argued that for successful operation 
in the quantum regime (i.e. when single atomic excitations are 
stored) the n pulses would have to be precise to of order 1 /TV, 
where TV is the number of atoms. Typically, solid-state ensem- 
bles contain of order 10 7 to 10 9 atoms, such a level of preci- 
sion would thus be completely out of reach. The arguments of 
Ref. i24ll were criticized in Ref. 01, but to our knowledge the 
question has not been fully resolved until now. 

In the present paper we study this problem in detail. The 
argument of Ref. Il24ll was based on the fact that imperfect 
ji pulses will lead to unwanted atomic excitations, which will 
cause background emissions. However, we will see that these 
emissions are non-collective and hence non-directional. As a 
consequence, good memory operation is achievable with real- 
istic n pulses. 

There are several different ensemble-based quantum mem- 
ory protocols. For definiteness, in the following we focus 
on the well-known Duan-Lukin-Cirac-Zoller (DLCZ) proto- 
col yfl. However, with small modifications our results ap- 
ply to many other protocols, including storage based on elec- 
tromagnetically induced transparency [6], off -resonant Raman 
transitions Jg], controlled reversible inhomogeneous broaden- 
ing |0], and atomic frequency combs lfl3ll . 

In the DLCZ quantum memory protocol, as shown in 
FigQJa), an off-resonant write pulse undergoes Raman scat- 
tering, leading to the creation of a single photon and a single 
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FIG. 1: Basic level scheme in the Duan-Lukin-Cirac-Zoller proto- 
col [3fl- ( a ) The far detuned write pulse scatters a write photon and 
creates a single collective atomic excitation (spin-wave) in the state 
s. (b) Applying a n pulse on the g — s transition in the middle of 
the process interchanges the roles of g and s, leading to rephasing at 
the end of storage period, (c) Shining the read pulse transforms the 
single collective atomic excitation in g into a read photon. 



collective excitation in the state s. This collective excitation 
dephases due to inhomogeneous broadening of the g — s tran- 
sition, but the application of a n pulse in the middle of the 
storage time, see FigQJb), can prepare a rephased atomic col- 
lective excitation at the time of retrieval. A read pulse can now 
be applied, which leads to the directional emission of the read 
photon, see Fig. [Tf c )- 

We study the effects of spin-echo related imperfections in 
the described protocol. We begin with a semi-classical treat- 
ment for uniform errors in the n pulses in section II. Then 
we give a fully quantum-mechanical treatment for the (most 
relevant) case of small n pulse errors in section III, and we 
show that its results agree with the semi-classical approach. 
We treat the case of non-uniform n pulse errors in appendix 
A. In section IV we consider the effects of imperfect optical 
pumping, i.e. an imperfect initial state. In section V we use 
the semi-classical approach to discuss the application of mul- 
tiple n pulses (bang-bang control). In section VI we give our 
conclusions. 



II. SEMI-CLASSICAL APPROACH 

We start with a semi-classical approach. Here the collec- 
tive atomic state is treated as a tensor product of single-atom 
states. The single-excitation component of this tensor prod- 
uct corresponds to the true quantum state, which is why the 
two approaches give equivalent results in the single-photon 
regime, cf. section III below. 

Consider an ensemble of A type three-level atoms with two 
slightly split ground states g and s, and an excited state e. 
At the beginning, all of the atoms are ideally pumped into 
the ground state g. Applying the write pulse that scatters a 
(Stokes) photon, transforms the state of the k' h atom at the 
position X/t into, 



l<A W (fo)> = \g) ~ & 



s), 



(1) 



where £ represents the contribution of each atom to the col- 
lective excitation, and Aki = k w - k s , where k w is the ^-vector 



of the write pulse and k s that of the scattered (Stokes) photon. 
For single-photon storage using an ensemble with N atoms, 
Ng 2 = 1. The total atomic state, as the product of single 
atomic states, implies that in the limit of £ — > the semi- 
classical atomic state tends to quantum mechanical state, see 
below. For weak light storage this limit is equivalent to limit 
of large number of atoms, N » 1 . 

Now, let us consider the effect of inhomogeneous broad- 
ening and the n pulse. In the semi-classical picture one can 
easily use unitary evolution of a single atom in absence of the 
electromagnetic field to represent the dephasing (rephasing) 
before (after) applying the n pulse. The propagator is given 
by, 



U Ak (t f ,td 



1 

e -' A *('/-'i) 



(2) 



where is the detuning from the central transition for the 
k ,h atom. The At have an inhomogeneous distribution with 
a width F. Then, the state of an atom after the time interval 
ti = fi - t is \iff (k \ti)) = U Ak (t u t )\^ k) (t )). An efficient re- 
trieval is impossible for long storage times such that Ft > 1, 
since the readout amplitude is governed by the average atomic 
polarization, which is greatly reduced by the dephasing that 
takes place due to the temporal phase factors e~' A '( f /~';) that 
vary from atom to atom. Applying the rf n pulse, which is 
tuned to the central frequency of the g - s transition, can bring 
this random phase into a negligible global phase at a certain 
time. In order to represent the n pulse, let us recall the ex- 
pression of the propagator of two levels of the atom under a 
pulsed excitation with the Rabi frequency £2, in the rotating 
wave approximation, 



U"(T) = 

9; = SliT, 



cos(6»,/2) -zsin(0;/2) 
-i'sin(0;/2) cos(0 ; /2) 



(3) 



where T is the temporal duration of the pulse. 

The final state after applying a ji pulse at t — t\ and 
waiting the time interval T2 = t<i — t\ is It^'fe)) = 
U At (t 2 , ti)U (T)U At (t u f )|iA w (fo)), where 6 = n+e and6» = n 
represents a perfect n pulse and e is the error. One can now 
retrieve the read photon by applying the read pulse. Hence, 
we need to consider the spatial phase that comes from the 
read pulse to find the direction of the read photon. The spatial 
phase dependence due to the rf pulse can be ignored, because 
krf.Xw <K 1 for realistic sample dimensions. Considering this 
fact, the final state after applying the read pulse is given by, 

\tj/®) = e £ ' lt (cos(6/2)-^sm(0/2)e- iAkT 'e iAl ' jlt )\e) (4) 



je- iAtT2 (sin(6»/2) + £T iA * T| cos(6»/2)e ,Al ' ft )|s), 



where k r is the fc-vector of the read pulse. 

By transferring the population of the state g into an excited 
state, the spin coherence transforms into an optical coherence, 
which leads to emission of the optical echo. The following 
shows how the atomic polarization would serve as the source 
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FIG. 2: A schematic representation of the collective enhancement 
at a certain direction that is given by the phase-matching condition. 
Non-directional re-emission is suppressed by a factor 1 IN, where N 
is the number of atoms. 



of the echo signal, 



*echo 



/() 



\Pf\ 2 



(5) 



P f = ^(e \ i/ffWp | s), 



where, y. is the electric dipole moment and Iq is the radia- 
tion intensity of one isolated atom. It can be seen from eq. 
© that terms with atom-dependent temporal phases appear 
in the atomic polarization in the eq. (0. Since the emis- 
sion is governed by the average over single atomic polar- 
izations, only those terms for which these phases are can- 
celed contribute significantly to the echo intensity. Analyz- 
ing different terms in (e \ i/r?)(if/j) \ s) one finds that the 

term ^sin 2 (6'/2)e' Al(T ^ T ' ) e , ' (Al ' 1+ ^ ) ;?1 is the only one for which 
one can exclude the atom-dependent temporal phase, which 
is called rephasing. This can take place under the condition 
T\ = T2, which implies that the n pulse has to be applied at the 
middle of the process. 

The dipole moment of each single atom serves as a source 
of radiation. In the far field approximation and under the 
rephasing condition (ti = T2), one can show that the am- 
plitude of the readout from the whole ensemble in the di- 
rection k ro is proportional to Y/J =i e «(A*i+Afe)-*j s i n 2 g/2, where 

A&2 = k r — k ro . Hence, for Ak\ + A&2 = the readout intensity 
is proportional to N 2 . This corresponds to constructive inter- 
ference of the re-emission from all of the atoms at a certain 
direction k ro - k w + k r - k s . However, even by applying an 
ideal n pulse the radiation at other directions is not zero. The 
background re-emission (non-directional) intensity is propor- 
tional to =1 e^i+AkM-l*)) = N, where A\ + Ak 2 * 0. 
Accordingly, the intensity ratio of the collectively enhanced 
re-emission (directional) and the randomly distributed back- 
ground re-emission (non-directional) is N, see Fig. ©. For 
example, for counter-propagation of the write and read pulses 
the phase-matching condition gives k ro = —k s . 

However, the error in the n pulse prevents achieving the 
highest possible the echo amplitude and presents a source of 
noise. 

Efficiency reduction. By considering these points one can 
calculate the atomic polarization, Pf - -iN^ji cos 2 (e/2), be- 



cause — k ± e. Consequently, the intensity of the echo is, 

I ec ho = IoN 2 fcos 4 (e/2), (6) 



where Iq has the same definition as in eq. (|3}. Eventually, the 
total efficiency of any process depends on the optical depth, 
performance of the experimental facilities and other theoreti- 
cal and experimental details which are related to the protocol 
and the experimental setup. However, the present result al- 
lows to find the efficiency reduction due to the error in the n 
pulse that is important in our analysis. 

Noise-to-signal ratio. Due to error in the rephasing pulse 
(e + 0) there is a noise in the re-emission. It is important to 
analyze the intensity of the noise, because for high memory 
fidelity the noise to echo ratio has to be small. 

A single n pulse interchanges the states of the atoms from 
s to g and vice versa. However, an error in the n pulse pro- 
duces population in the ground state g even without apply- 
ing write pulse at the beginning. This would lead to fluores- 
cent radiation. Hence, it can be distinguished by the non- 
zero terms in ') for £ = 0, which means even with no 
write pulse the error by itself can produce population in g. 
Quantitatively, the fluorescent radiation can be specified by 

I (g I *A/ ) ^°> l 2 = cos 2 (6/2) which is the term in \\Jt f ) that 
does not originate from the read-in process (^ = 0), but gives 
a non-zero projection on \g). Therefore, the intensity of the 
noise that originates from error in the n pulse can be given by, 



(7) 



k=\ 



/ JVsin 2 (e/2). 



Obviously, the fluorescent radiation as the source of the noise 
has the spatial dependence of a single photon radiation. Con- 
sequently, the semi-classical picture yields an equally dis- 
tributed noise (non-directional) that comes from single atom 
radiation of ensemble of the atoms, because of error in the 
rephasing pulse. 

In order to analyze the effect of an imperfect 71 pulse on the 
fidelity, it is important to study the noise-to-signal ratio. The 
following presents the noise-to-signal ratio, 



I mise sin (e/2) 



*echo 



cos 4 (e/2)' 



(8) 



keeping in mind that A^^ 2 = 1 for single-photon storage. The 
higher the noise-to-signal ratio, the less fidelity we have. In 
112311 . 1% variation has been realized in the intensity of the rf 
pulse that causes the same error in the pulse area. Such error 
gives quite low noise-to-signal ratio of 0.25 x 1CT 4 and only 
wastes 0.1% of the efficiency. The semi-classical calculation 
suggests that the typical 1 % error in the n pulse which is far 
beyond the 1 /N precision, does not impose a major constraint 
on the efficiency and fidelity. In addition to the semi-classical 
calculations, it is interesting to reconsider the problem using 
the quantum mechanical description. In the next section, we 
perform fully quantum mechanical investigation for the same 
question. 
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III. QUANTUM MECHANICAL TREATMENT 

In this section, we perform fully quantum mechanical anal- 
ysis on a global error in the n pulse. We consider an ideal 
initial state where all the atoms have been pumped into the 
ground state. 



A. Ideal Protocol 

As we discussed before, the large detuning of the write laser 
from g—e leads to scattering of a photon and creating a collec- 
tive atomic excitation. The theory of light-atom interaction is 
well established to describe interaction between the field and 
the atomic dipole moment. In general, one can represent the 
interaction Hamiltonian as, 



N r 



dka(k)e' kx 'd J pv + h.c, 



(9) 



where G = (p\p-j.e\v) and & J pv = |p> 7; (v|, where p and v 

denote atomic levels g, e and s. This Hamiltonian shows the 
interaction between a field and the dipole moment of the f h 
atom, fij. For simplicity, we assume the dipole-field coupling 
is identical for all the atoms and we have not considered the 
transverse profile. One can extract the interaction Hamilto- 
nian for the write laser and the scattered photon from eq. (0, 
see FIGQJa). By combining the interaction Hamiltonians and 
considering adiabatic elimination of the excited, <?, one can 
describe the read-in part of the process as follows, 



H 



■ff 



N 

z 



G' j dk s al(k s )e i{{ " 



u sg 



■ h.c. 



(10) 



where G' = (s\flj.? s \e)(e\fij.e w \g} J 0e».(r) and Q h .(t) = 

(e\p.j.e„.\g)s w (T) is the Rabi frequency of the classical write 
field. 

The unitary evolution under this effective interaction 
Hamiltonian describes the creation of a Stokes photon via Ra- 
man scattering, which is accompanied by the creation of a 
collective atomic state. The collective state is 



Wo)> = 



1 



N 



Vat ^ 



\g.-s {k) .:g), 



(ii) 



i<=\ 



where \g...s^ k \..g) shows all atoms in the ground state and the 
k' h atom in the spin state s and Aki = k w -k s . Here, we neglect 
the terms with more than one excitation, because we are inter- 
ested in single photon storage. The interaction Hamiltonian 
governs the evolution of the atomic and the photonic state of 
the system. Our analysis follows the evolution of the atomic 
state. However, later in this section, we will refer to the above 
discussion in order to find the intensity of the readout based 
on the norm of the final atomic state. 

One can find the state \if/(to)) in eq. ( fTTT > by apply- 
ing the Schwinger bosonic creation operator or equally 



weighted superposition of <t+ =| s)(g | operators, J + (Aki) = 

2^ = i e iMl 'erf <g> 1. The inhomogeneous spin broadening im- 
plies that each atom has a slightly different energy than the 
other atom's spin level that indicates any atom will evolve 
based on its small energy detuning specified by A„. After the 
time interval t\ — t\ - to the state |i^(fo)) will be evolved to 
\i//(ti)) = A=J^ =x e ihkTi e i&lSk \g...s ( - k \..g). The effect of de- 
phasing as a result of inhomogeneous spin broadening can 
be described as an atom-dependent phase accumulation of the 
single collective atomic excitation. 

The Following operator is appropriate to mathematical 
modeling of the dephasing after the time t, 



e iat = , 



"k=l 



(12) 



where £2 = Y*k=i Ak\s)k{s\ ® 1, see [26], where the operator 1 
shows the identity operator that acts on the rest of the atoms. 
The operator in Eq. (TT2l can be used to show the rephasing 
after applying the n pulse. Hence, one can represent the de- 
phased state as |0fo)> = -j=e ihT 'J + (Ak 1 )\gg...g). 

At first let us recall how an ideal n pulse treat the atoms, 
which it swaps g to s and vice versa. One can describe that as 

e m/2J x = e m/2^ «^.i > where ^pj, = ^A^W see ^ 

After applying the ti pulse, the operator e' n ' leads to a 
rephasing of the collective excitation because for each term 
in the collective state the previously non-excited atoms now 
acquire phases whereas the previously excited atom doesn't. 
Finally, the stored pulse will be retrieved by applying the 
read field, leading to the emission of the readout photon, 
see FIGOIc). The readout part of the process can be repre- 
sented by the operator J+(Aki) = 2f=i e'^ kl ^ k cr ( + ® 1, that is 
based on the same discussion as for the read-in. Eventually, 
the whole process comprised of creation of collective exci- 
tation, dephasing, n pulse, rephasing and the read pulse that 
leads to the final state can be described as \ifff(Ak\, A&2)) = 
-^J + (Alc 2 )e^e^ 2J ^J + (Alc 1 )\gg...g). 

Using Eq. ( TTOb and the analogous Hamiltonian for the read- 
out process it is straightforward to include the quantum states 
of the light field for the Stokes and readout photons into the 
description. One sees that the emission amplitude for the 
readout photon can be obtained directly from the norm of the 
atomic state |0y(Afci, Afei)). Note that there is no preferred di- 
rection for the Stokes photon emission |27]. As before, we 
are not really interested in the absolute emission probability, 
but in how the probability varies as a function of the direction 
of emission for the read photon, and under the influence of 
errors in the control pulses. The collective enhancement hap- 
pens again under the phase-matching condition Aki +Ak2 = 0. 

In the ideal case, considering these conditions one can eas- 
ily derive that {>J/f(A\, A\)\ift f (Ak u Ak 2 )} = N. The later 
result is based on considering an ideal n pulse in the calcula- 
tions and it takes place under the phase-matching condition, 
A£i + A&2 = 0. As we studied in the semi-classical approach, 
the re-emission in other directions {Ak\ + Ak 2 + 0) is negligi- 
ble for the case with large number of atoms in the ensemble. 
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Our discussion here was focused on the case of the DLCZ 
protocol, but the evolution of the atomic state is identical in 
the other quantum memory protocols mentioned in the intro- 
duction. The only difference is that in those protocols the cre- 
ation of the initial atomic excitation is accompanied by the 
absorption of a single photon rather than its emission. 



B. Global error in the n pulse 

We now study the effect of imperfect n pulse as one of the 
sources of error in the spin echo memories. The most gen- 
eral error in the n pulse can be considered as a small rotation 
€k around a random direction hk for each atom that can be 
represented by the operator e^*=' "*/ 2o4 "i ai| »i xhis operator 
can show the effect of the inhomogeneity in the intensity of 
the n pulse across the sample ifzll . For simplicity, we first 
consider a global error that can be interpreted as lack of accu- 
racy in pulse shaping. The most general error that affects each 
atom differently is studied in the appendix. The correspond- 
ing operator for such error which is equally distributed over 
all atoms is e iel2]A = e ie/2 Zli <?'"-"®i. 

Essentially, the question of studying the effect of the error 
in the n pulse is reduced to the problem of analyzing norm of 
the state 

|^/(A^i, A^ 2 , e)> = ^/ + (Al 2 )Av >/2 V e/2y " (13) 
e ihT 'J + (Ak l )\gg...g}. 

In order to facilitate analyzing the norm of the state 
\if/f(Ak\, A&2, e)) we now simplify the final state. Since we 
are only interested in its norm, any unitary operator can be 
used to simplify the state. By applying the unitary operator 
e -m/2J x f rom the left and adding the identity e m ^ 2J 'e^ m ^ 2J -\ one 
can represent the final state up to a unitary and a global phase 
as, 



1 



J-(A\)e- in ^-e ie/2J!> e inT >J + (A\)\gg. 



(14) 



where J_(Ak2) 



,iA 2 l k(J ik) 



1. Under the condition 



t — T\ — T2 and by conducting some algebra one can derive 
that e -^T e iei2j.h e iCir _ e zL'V2#* ) .»'*>®i ) which shows rota . 

tion in a new direction h' (k) = (ji^cosAjt + n v sin A^t).* + 
(-n x sinAAT + n y cosA J tT)y + n z z. Finally, by applying the 
unitary e^*= ] -< e t/2<x m .«' ( "«i me fj na j state can ^ e simplified to 
OJ+\gg...g) where the operator O is given by, 



6 = -4= y e^iao*? + y8e iA *Vf + ye 2iA *Vf ) ® 1(15) 



k=] 



where a = cos 2 (e/2) + 2/ sin (e/2) cos (e/2) - n 2 sin 2 (e/2), 
/3 = -i sin (e/2) cos (e/2)(n x - in y ) + n z sin 2 (e/2)(n x - in y ) and 
y = sin 2 (e/2)(n x - in y ) 2 . This simplification allows us to dis- 
tinguish three different terms in the final state based on the 
effect of OJ+ on the \gg...g). The first term corresponds to 



cr^<T ( ^\gg...g) which gives 6jk\gg—g). The other combina- 
tion erf V!^ in OJ+ yields (-l) s > k \g..s u) ..g).The later gives N 2 
terms that contains one excitation. Finally, the erfcr*/' leads 
to N(N - 1) terms with two excitations. One can benefit from 
these terms to represent the final state as 



a|i/ri)+/3|lA2> + rl<A3>, 



(16) 



where |^> = 4- Z? =1 e* A ^ 1 i\g..g), Mi) = 



^I^ M=l {-l) 5lk e i ^ T e^ ll ie i& * i <'\g..s <J) . 



and |(/r 3 > 



-J^'Lj,k=ij*k e2iAkTe ie I g»s U) ..s®..g). Obviously, 
because of different number of excitations, these terms are 
perpendicular. Thus, one can study the norm of the final state, 
easily. 

As it can be seen from the eq. ( TTBT l. because of Sjk, the first 
term corresponds to the directional emission of the readout 
photon. As we discussed in the previous section, see FIG. 
(fJJ, the readout is strongly peaked around the direction for 
which all of the single-atom re-emissions can constructively 
interfere with each other. The readout intensity at the other 
directions (the non-directional background) is suppressed by 
the ratio 1 /N that is a quite small number for a typical N ~ 
10 7 - 10 9 . Hence, the a term determines the intensity of the 
echo. In other words, shining the read laser with k r = -k w 
result in the emission of a readout photon that its intensity is 
peaked around k ro = -k s . The non-directional noise in the 
re-emission can be attributed to the terms that correspond to 
/3 and y coefficients, because taking average over position of 
the atoms leads to randomly distributed re-emission from each 
single atom. It can be shown that the norm of the state in Eq. 
©is 

$ (17) 

N N . N 

7=1 *=1 N j,k=l,j*k 

Efficiency reduction. In order to study the worst case sce- 
nario, upper bound of the noise strength and lower bound 
of the echo intensity have to be studied separately. Con- 
sidering the phase-matching condition it can be shown that 
hcho N\a\ 2 . For small errors, e <K 1 and keeping terms up 
to 0(e 3 ), it can be shown that the following lower bound can 
be achieved for n, = 0, 



|a| 2 <(l-2(e/2) 2 ). 



(18) 



As in the ideal case the signal intensity is proportional to 1, 
one can analyze the efficiency reduction. So from quantum 
mechanical point of view the efficiency reduction factor is 
given by 1 - 2(e/2) 2 that gives the same results for the typ- 
ical 1 % error as the semi-classical calculation. 

Noise. With the aim of studying the in- 

tensity of the noise that is proportional to the 

^l/?l 2 2£, I ^(-ly^e^e^? + 1 X% =um \y\ 2 
in the limit of the small errors, e <k 1, we need to discuss the 



6 



effect of e kT . For long enough times that t is comparable 
with i, where F is the inhomogeneous linewidth, the e' A ' T 
becomes a completely random phase. It can be shown that for 



random phases 0*, | e*Q| — I 2*,/ Qc ; *e 



'(01- -0/) I 



because (<? !( * = 6ki- Hence we can treat the term that cor- 
responds to p as j;\P\ 2 lfj,k=i K-l^e^i^g^^p = N\p\ 2 . 
Taking these considerations into account shows that the 
intensity of the noise is proportional to 

1 N N 

-\ p | 2 Y| V(_l)^ e iA t r e i(Al 1 )4 e i(Al 2 ).Z t |2 (19) 

7=1 *=1 
#(#-1). |2 



+ 



N 



-Irr 



< JV(e/2) 2 max(n 2 (n 2 x + n 2 ,)(e/2) 2 + n 2 + n 2 ) 



N 2 - N 



(e/2) 4 max|(n,-m y ) 2 | 2 



iV 

~ Me/2) 2 + <9(e 4 ). 

So the choice of uniformly directed error with n, = gives the 
upper bound for the noise. Obviously, the second term in the 
noise intensity is proportional to e 4 that is negligible for the 
small errors. Then it leads to I no j se °c N(e/2) 2 . These results 
allow us to find the upper bound for the noise-to-signal ratio. 

Let us recall the results for the noise and the echo from the 
semi-classical approach in eqs. ( 1617b . Obviously, the Taylor 
expansion of the results obtained in semi-classical treatment 
and eliminating the terms 0(e 4 ) and higher, demonstrates the 
agreement between the results of the both approaches in the 
limit of small errors (e <k 1). 

In the quantum mechanical approach the noise is propor- 
tional to the term given in eq. ( fT9l ). The amplitude of the 
noise shows a correspondence with the fluorescent radiation in 
the semi-classical approach. Indeed the spatial phase depen- 
dence implies that the direction of the emission of the noise 
from each atom varies from one to another, leading to a non- 
directional noise. Thanks to collective enhancement, the non- 
directional noise will not swamp the echo signal, for realistic 
control pulse accuracy l23ll . 



IV. IMPERFECT INITIAL STATE AND n PULSE 

An imperfection in the n pulse is not the only source of inef- 
ficiency in the spin echo memories. In our calculations so far 
we have assumed an ideal situation where all the atoms are ini- 
tially in the ground state. Our quantum mechanical approach 
allows us to study the effect of an imperfect initial state with n 
atoms excited to the state s. This can happen in experiments as 
a result of an imperfect optical pumping in the initialization of 
the atomic ensemble. Without loss of generality, we can con- 
sider the initial state \g..gs..s) that has n sorted excited atoms 
instead of randomly positioned excited atoms. Applying the 
operator OJ+ gives the final state. We expect the directional 
echo as result of applying Yj N jk=\ o"_^cr+ on the initial state 
\g..gs..s). In contrast to the perfect initial case, this will lead 
to N - n terms with n excitations correspond to \g..gs..s), and 



also (N - n)n terms connected with \g..s^\.gs..g^ k \.s) which 
has n excitations. By conducting some algebra, one can eas- 
ily show that only the first case gives Sjk which leads to di- 
rectional (collectively enhanced) re-emission under the phase- 
matching condition. This implies that by considering imper- 
fection in both the initial state and the n pulse, intensity of the 
echo is proportional to 



lecho O*- 



(N - n y 

N 



(l-2(e/2) 2 ). 



(20) 



The imperfection in the initial state reduces the intensity of the 
echo and also introduces a new source of the non-directional 
noise. Previously, we analyzed the noise that corresponds to 
\if/2) and It/^) by finding the upper bound for /3 and y. Now, the 
(N-ri)n terms of the form \g..s ( j\.gs..g( k \.s) also contribute to 
the noise. Consequently, as the upper bound, for the term that 
corresponds to a can be achieved for n z = 1, one can obtain 
the following upper bound for the noise-to-signal ratio r, 



r < 



(N - n)n + N 2 (e/2) 2 
(/V-n) 2 (l -2(e/2) 2 )' 



(21) 



Investigating this equation it can be seen that the two different 
sources of error compete in producing noise. For the case with 
small fraction of excitations in the initial state, ^ <K 1, such 
that e 2 > y then the error in n pulse is the dominant term in 
the noise. Otherwise, in case of ^ > e 2 the imperfection in 
the initial state plays an important role in increasing the upper 
bound of the noise-to-signal ratio. For instance, 2% of the 
atoms not in the ground state, which is a typical value fiUl . 
will cause a 2% error. The error in the n pulse would have to 
be as large as e = 0.4 in order to be comparable in importance. 



V. SEMI-CLASSICAL TREATMENT OF MULTIPLE n 
PULSES 

So far we have studied the application of a single n pulse. 
Much longer storage times are achievable by applying se- 
quences of n pulses. This is also known as bang-bang control. 
We analyze this case using the semi-classical approach. We 
have seen before that semi-classical and quantum treatment 
lead to identical conclusions for small errors. The intensity of 
the signal after applying m pairs of n pulses is given by, 



lecho 



hN 2 f(l 



(2m 2 -m + 1)^12). 



(22) 



The intensity that is associated to the noise also can be approx- 
imated by l no ise ~ IoNm 2 £ 2 /2. Thus, after applying sequence 
of the rephasing pulses, the noise-to-signal ratio reads 



,2,2 



e l l2, 



(23) 



for me <K 1 . For instance, with typical 1 % error in the n pulse 
one can benefit from 30 pairs of rephasing pulses, while still 
achieving an efficiency factor of 91% and a noise-to-signal 
ratio of 0.04. If the pulses are 4 msec apart as in Ref. 10], 
30 pairs of pulses correspond to a storage time of 240 msec. 
Ref. rt23ll managed to increase the T2 time to 1 sec, by proper 
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alignment of the magnetic field, that is much longer than 80 
msec in Ref. [!7]. This field configuration would allow one to 
put the rephasing pulses further apart. For example, one can 
consider the rephasing pulses 40 msec apart, giving a 2.4 sec 
storage time with an efficiency factor of 91% and a noise-to- 
signal ratio of 0.04. 



VI. DISCUSSION AND CONCLUSION 

We have shown that realistic imperfections in the n pulses 
only have small effects on the retrieval efficiency of quantum 
memories, and that the corresponding noise is also accept- 
able. The latter fact is due to collective interference. While 
the emission of the read photon is strongly enhanced in the di- 
rection given by the phase matching condition, the noise due 
to pulse imperfections is non-directional. As a consequence, 
even the use of great numbers of n pulses (bang-bang con- 
trol) is realistic at the quantum level. This settles the question 
raised by Ref. B2411 . 

We have also studied errors due to imperfect optical pump- 
ing, i.e. imperfect initial state preparation. We find that these 
errors are likely to dominate over n pulse errors in many cir- 
cumstances, but that they can also be kept at an acceptable 
level. 

Our discussion was phrased in terms of the DLCZ proto- 
col JH for concreteness. However, the same results apply to 
other ensemble-based quantum memory protocols including 
those based on electromagnetically induced transparency [6], 
off-resonant Raman transitions |8], controlled reversible inho- 
mogeneous broadening Jgt], and atomic frequency combs lfl3ll . 

In conclusion, the prospects for the use of spin-echo tech- 
niques in light-matter interfaces at the quantum level are very 
good. We hope that our work will further encourage experi- 
mental work in this direction. 
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Appendix A: Spatial inhomogeneity in the rf pulse 

As we discussed, the main error in rf pulses is the variance 
in the intensity of the rf pulse that leads to inhomogeneity 
of the rf pulse across the sample II23I1 . Consequently, the er- 
ror would be different for the atoms at the different positions. 
Here we extend our analysis to study such errors. 

We study the spatial inhomogeneity in the rf pulse by con- 
sidering ' e </ 2cr A ®i as me operator which represents the 
error. Therefore, based on eq. ( fT~3l > the final state is, 



\tfr f (Ak u Ak 2 ,e)) = — — /+ (M 2 )e'" T2 e in/2J * 
V7V 



(Al) 



We can follow the same approach as we used in the paper 
to simplify the final state that eases the rest of the calculation. 
Therefore, the final state |0-/(AAri, AA:2, £)) up to a unitary op- 
erator and a global phase is 



1 



J-(Ak 2 )e- li 1T2 e 



J+(Ak] 



(A2) 



Considering the n pulse at the middle of the process, r = 
Ti = T2, and by conducting some algebra one can derive that 
e -i&r e Xl, i^.^i e inr = g zL which shows rota . 

tion in a new direction h' {k) = (n® cos A^t + riy sin A k T)x + 
(— sin AkT + nf 1 cos A k r)y + nf'z. Finally, by applying the 
unitary e^JLi l &°- m -n' {t) ®i me fj na j state can ^ e simplified to 



\^ f (A\ u Ak 2 ,e))=OJ + \gg...g), (A3) 



where the operator O, which acts on k' h atom is given by, 



1 N -, -, 

= — V + fF>e** T o? (A4) 

+ y^e^o-f) ® 1, 



and oP * 



cos 2 (6^/2) + 2/sin(e / t/2)cos(e / ( : /2) - 



nf sin 2 (e k /2), p {k) = -ism(e k /2)cos(e k /2)(nf - inf) + 
nf sin 2 (e k /2)(nf - inf) and y« = sin 2 (e t /2)(nf - inff. 
Finally, the state can be simplified to 



1 N 

^ f (Ak u A\,e)) = — J]a^e* Ak ^ 1 J\g..g)(A5) 



+ _!_ (-l) s ' t B (k) e i ^ kT e m ' )Xl e i ' 



j,k=\ 



Zf =1 ie t /2o'<«.n (t )®l iOr 



e* lT, J + {Ak\ 



1 N 

— ( y - 



?.J j \.s ik) ..g). 



Fortunately, having different number of excitations in the final 
state simplifies the calculation of the norm. 

As we discussed, the first term contributes in directional 
emission of the readout photon. Hence, the a^s play role in 
finding the intensity of the echo. In other words, shining the 
read laser with k r = —k w result in the emission of a readout 
photon that is peaked around k ro = -k s . The non-directional 
noise in the re-emission can be attributed to the terms corre- 
sponding to /3 (k) and y ik) coefficients, because average over 
the position of the atoms leads to randomly distributed re- 
emission. It can be shown that the following gives norm of 
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the \4ij) in eq. ( lA5l ). 



1 N 



7=1 

7V TV 



0')^i(Ali+A^ 2 ).^ |2 



(A6) 



+ i V | ^(_ 1 )^ e ''A i T^) e /(Al 1 ).f Je ;(Al 2 ).f t ,2 
7=1 t=l 
1 W 

+ 1 y 

A/ Z-l 



Jk) ,2 



In order to analyze the efficiency reduction the lower bound 
of echo intensity have to be studied. Considering the phase- 
matching condition it can be shown that I ec i,„ oc jj\ Y^_j a ,(; '| 2 - 
For small errors, ej <s 1 and keeping terms to 0(e 3 ), it can 
be shown that the following lower bound can be achieved for 



n? = 



1 N 

-| J] aU) \ 2 * N(l ~ 2(w2/) 2 ), (A7) 

7=1 

where e max is the largest error and N is number of 
the atoms. With the aim of studying the inten- 
sity of the noise that is proportional to the ^ Tlj=\ I 

2f =1 (-l)^e iA ^e« A ')^e^)^ | 2 +i^ M#i | r « | 2 



for the small errors, e, <k 1, one needs to consider the e' AkT 
as a random phase, that takes place for long enough times that 
r is comparable with the p. It can be shown that for ran- 
dom <i, k , iEt«*jS»f = \Zk,iP {k ¥' r ^- M \ 2 = Zk\P (k) \ 2 , 

because (e^'W) = 5 kl . Hence, one can conclude that 
I 2f =1 (-l) 5 ' 1 e , ' A4T /3 w e i(All)f 'e i(Al2)f '| 2 can be approximated as 
ZL |(-l)^W e «^^ e ; (AlM|2 = 2f =1 |(-1)^/JW|2. By 
taking these considerations into account and keeping terms up 
to 0{e A ) give the noise intensity upper bound as 



(A8) 



< JVmaxCnf (nf +nf ) )(e k /2? 

+ («f 2 + »f)(e t /2) 2 ) 

+ —Z-^maxCC^KfiW -mf) 2 | 2 ) 

= jV(e fflfl ,/2) 2 +0(e 4 ), 



which implies that a uniformly directed error with n, = gives 
the upper bound for the noise. Obviously, the second term in 
the noise intensity is proportional to e 4 that is negligible for 
the small errors. Then it leads to I„ ojse < N(e max /2) 2 . This 
shows that the fully quantum mechanical treatment for an in- 
homogeneous error is in good agreement with the results for 
global errors derived in the paper. 
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